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Abstract 

It is shown that in the case of the one-particle one-dimensional scattering prob- 
lem for a given time-independent potential, for each state of the whole quantum 
ensemble of identically prepared particles, there is an unique pair of (subensemble's) 
solutions to the Schrodinger equation, which, as we postulate, describe separately 
transmission and reflection: in the case of nonstationary states, for any instant of 
time, these functions are orthogonal and their sum describes the state of all particles; 
evolving with constant norms, one of them approaches at late times the transmitted 
wave packet and another approaches the reflected packet. Both for transmission and 
reflection, 1) well before and after the scattering event, the average kinetic energy 
of particles is the same, 2) the average starting point differs, in the general case, 
from that for all particles. It is shown that for reflection, in the case of symmetric 
potential barriers, the domain of the motion of particles is bounded by the midpoint 
of the barrier region. We define (exact and asymptotic) transmission and reflection 
times and show that the basic results of our formalism can be, in principle, checked 
exp er iment ally. 



PACS numbers: 03.65.Ca,03.65.Xp 
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1 Introduction 



For a long time tunneling a particle through an one-dimensional time-independent poten- 
tial barrier was considered in quantum mechanics as a representative of well-understood 
phenomena. However, now it has been realized that this is not the case. The inherent 
to quantum theory standard wave-packet analysis (SWPA) [1, 2, 3, 4, 5] (see also [6]), in 
which the study of the temporal aspects of tunneling is reduced to following the centers of 
"mass" (CMs) of wave packets, does not provide a clear prescription both how to interpret 
properly the scattering of finite in x space wave packets and how to introduce charac- 
teristic times for a tunneling particle. All these questions constitute the main content of 
the so-called tunneling time problem (TTP) which have been of great interest for the last 
decades. 

As is known, the main peculiarity of the tunneling of finite wave packets is that the 
average particle's kinetic energy for the transmitted, reflected and incident wave packets 
is different. For example, in the case of tunneling through an opaque rectangular barrier, 
the average velocity of the transmitted particle is larger than that of the incident particle. 
It is evident that this fact needs a proper explanation. As was pointed out in [7, 8], 
it would be strange to interpret the above property of wave packets as the evidence of 
accelerating a particle (in the asymptotic regions) by the static potential barrier. Besides, 
in this case there is no causal link between the transmitted and incident wave packets 
(see [7, 8]). 

As regards wide (strictly speaking, infinite) in x space wave packets, the average kinetic 
energy of particles, before and after the interaction, is the same. But now the uncertainty 
in defining the CM's position and, consequently, corresponding asymptotic times is very 
large. As a result, the most of physicists considers the characteristic times introduced in 
the SWPA as quantities having no physical sense. The review [1] devoted to the TTP 
seems to be the last one in which the SWPA is considered in a positive context. 

Apart from the SWPA, in the same or different setting the tunneling problem, a 
variety of alternative approaches (see reviews [1, 8, 9, 10, 11, 12] and references therein) 
to advance different characteristic times for a tunneling particle have also been developed. 
Among the alternative conceptions, of interest are that of the dwell time [13, 14, 15, 16], 
that of the Larmor time [17, 18, 19, 20, 21] to give the way of measuring the dwell time, 
and the conception of the time of arrival which is based on introducing either a suitable 
time operator (see, for example, [22, 23, 24, 25, 26]) or the positive operator valued 
measure [11]. Besides, of interest are attempts to study the temporal aspects of tunneling 
in the frame of the Feynman, Bohmian and Wigner approaches to deal with the random 
trajectories of particles (see, for example, [27, 28, 29, 30, 31] and references therein). One 
should also mention the papers [32, 33, 34] where the TTP is studied in the framework of 
a nonstandard setting the scattering problem. But again, for a particle whose initial state 
is described by a Gaussian-like wave packet, none of the alternative approaches have not 
yet led to commonly accepted characteristic times (see the reviews [1, 8, 9, 10, 11, 12]). 

Note, unlike the SWPA, in all these approaches theoretical efforts have been aimed 
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at elaborating some rule of timing the motion of a quantum particle. As was said in 
[11], "... up to now the interest of theoreticians [to the TTP] has been motivated ... 
by a fundamental lacuna of quantum theory, namely the absence of a clear prescription 
to incorporate time observables into its formalism" . However, in our opinion, there is no 
necessity in such a prescription: time in quantum mechanics has the status of a parameter, 
and, hence, there is no room here for time observables (or, time operators). And, what 
is more important, the rule of timing the particle's motion have already been available in 
quantum theory, and this rule is dictated by the correspondence principle. 

By the analogue with classical mechanics where timing the particle's motion is reduced 
to the analysis of the function x(t) (x is the particle's position, t is time), in quantum 
mechanics characteristic times for a particle should be derived from studying the temporal 
dependence of the expectation (average) value of the position operator for a particle in a 
given state (or, what is equivalent, from studying the temporal behavior of the CM of the 
corresponding wave packet). Besides, from the analysis of the temporal dependence of the 
mean-square deviation for this operator (or, from the analysis of the temporal behavior 
of the corresponding leading and trailing edges of the wave packet) one can take into 
account uncertainty in the particle's position, and thereby evaluate the error of the above 
timing. 

However, one has to bear in mind the following. The above timing procedure suggests 
that the average value of the position operator has its primary physical sense, as the 
most probable position of a particle, at all stages of its motion. For a free particle 
whose state is described by a Gaussian-like wave packet, this requirement is fulfilled and, 
thus, no problem arises in timing its motion. An essentially different situation takes 
place in the case of a tunneling particle. Now, following the CM of the wave packet to 
describe the state of the whole ensemble of particles becomes meaningless at some stages 
of scattering. In particular, after the scattering event, when we deal, in fact, with two 
scattered (transmitted and reflected) wave packets, the averaging over the whole ensemble 
of particles has no physical sense. Of course, in this case there is a possibility to define the 
individual average positions of transmitted and reflected particles. However, in timing, 
such averaging suggests the separate description of both the subensembles at early times, 
what is widely accepted to be impossible in quantum mechanics. As was pointed out in 
[12], "... transmission and reflection are inextricably intertwined". 

So, quantum theory provides the needed rule of timing the motion of a particle, but 
it conflicts with the existing viewpoint that transmission and reflection are allegedly in- 
separable (namely this obstacle have remained to overcome in the SWPA). However, in 
our opinion, there is no basis for such a viewpoint. For none of the principles of quantum 
mechanics forbids such a separation. In reality, the main problem is that quantum theory, 
as it stands, does not provide the way of separating transmission and reflection. In our 
opinion, just the absence of the corresponding mathematical formalism is a fundamental 
lacuna in quantum theory. 

In this paper, in the framework of the conventional quantum mechanics, we show 
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that, at least in one dimension, transmission and reflection can be described separately. 
It is surprising that such a description needs no innovation. As it has turned out, the 
separation of transmission and reflection is provided by the intrinsic property of the 
Schrodinger equation, which have been overlooked before. Namely, we show that in 
the standard setting the one-dimensional scattering problem for a given potential the 
Schrodinger equation possesses, in addition to the solution to describe the state of the 
whole ensemble of particle, an unique pair of other solutions which, as we postulate, 
describe separately transmission and reflection. The basis for such a postulate is that, for 
any instant of time, the subensemble's nonstationary-state wave functions are mutually 
orthogonal and their sum describes the state of the whole ensemble of particles; one 
of them causally evolves into the transmitted wave packet, and another approaches at 
late times the reflected one. The main peculiarity of stationary-state wave functions for 
transmission and reflection is that, for a given energy of particle, there is a point in 
the barrier region where these everywhere continuous functions have discontinuous first 
spatial derivatives. Nevertheless, this point is not a sink or source of particles for each 
subensemble. So that the norms of the corresponding wave packets are constant in time. 

Note, at present there is a paradoxical situation. Although the tunneling phenomenon 
have been known for a long time, the properties of tunneling proper have remained, in 
fact, unstudied. The point is that the "full" wave function to describe the state of a 
particle in the one-dimensional scattering problem relates to all particles of the quantum 
ensemble, rather than to transmitted particles only. In this connection, we hope that 
the formalism presented here will be useful for a deeper understanding of the tunneling 
process and, in particular, Hartman effect [3] widely discussed in the literature (see, for 
example, [35]). 

The paper is organized as follows. In Section 2 we pose a complete one-dimensional 
scattering problem for a particle, and display explicitly shortcomings to arise in the SWPA 
in solving the TTP. In Section 3 we present a renewed wave-packet analysis in which 
transmission and reflection are treated separately. In Section 4 we define the average 
(exact and asymptotic) transmission and reflection times and consider, in details, the 
cases of rectangular barriers and ^-potentials. 

2 Setting the problem for a completed scattering 
2.1 Backgrounds 

Let us consider a particle tunneling through the time-independent potential barrier V(x) 
confined to the finite spatial interval [a,b] (a > 0); d — b — a is the barrier width. 
Let its in state, \P in (x), at t = be the normalized function ^f^ t {x). The function 
is proposed to belong to the set consisting from infinitely differentiable functions 
vanishing exponentially in the limit |x| — > oo. The Fourier-transform of such functions 
are known to belong to the set as well. In this case the position and momentum 
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operators both are well-defined. Without loss of generality we will suppose that 

< >= °. < *!>I*!S< >= "o > o, < i^K}, >= n 

here Iq is the wave-packet's half-width at t — (Zo << a); £ and p are the operators of 
the particle's position and momentum, respectively. 

An important restriction should be imposed also on the rate of spreading the incident 
wave packet. Namely, we will suppose that the average velocity hko/m is large enough, so 
that in the incident wave packet its parts lying behind the CM within the wave-packet's 
half-width move toward the barrier. 

As is known, the formal solution to the temporal one-dimensional Schrodinger equation 
(OSE) of the problem can be written as e~ lHt / n ty i n (x) . To solve explicitly this equation, 
we will use here the variant (see [36]) of the well-known transfer matrix method [37] that 
allows one to calculate the tunneling parameters, as well as to connect the amplitudes 
of the outgoing and corresponding incoming waves, for any system of potential barriers. 
The state of a particle with the wave-number k can be written in the form 



* full (x; k) = A in (k)e ikx + B out (k)e- ikx , (1) 

for x < a, and 



* fvll (x; k) = A out (k)e tkx + B m {k)e^ k \ (2) 

for x > b, where A in (k) should be found from the initial condition; B in (k) = 0. The 
coefficients entering this solution are connected by the transfer matrix Y: 



which can be expressed in terms of the real tunneling parameters T, J and F, 



q = , exp [i{kd — J(k))] ; p = 
sjTik) 

T(k) (the real transmission coefficient) and J(k) (phase) are even and odd functions of 
k, respectively; F(—k) = n — F(k); R(k) = 1 — T(k); s = a + b. Note that the functions 
T(k), J(k) and F(k) contain all needed information about the influence of the potential 
barrier on a particle. We will suppose that the tunneling parameters have already been 
known explicitly. To find them, one can use the recurrence relations obtained in [36] just 
for these real parameters. 



R(k) 
Tik) 



exp 



7T 



+ F(k) - ks 



(4) 
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As is known, solving the TTP is reduced in the SWPA to timing a particle beyond 
the scattering region where the exact solution of the OSE approaches the corresponding 
in or out asymptote [38]. Thus, definitions of characteristic times in this approach can be 
done in terms of the in and out asymptotes of the tunneling problem. 

Note that in asymptote in the one-dimensional scattering problem represents an one- 
packet object to converge, well before the scattering event, with the incident wave packet, 
while out asymptote represents the superposition of two non-overlapped wave packets to 
converge, at t — > oo, with the transmitted and reflected ones. It is easy to show that 
in the problem at hand in asymptote, ty in (x,t), and out asymptote, ty out (x, t), can be 
written as follows 

* in (x,t) = -= \ f in (k,t)e lkx dk, f m (k,t) = A in (k)eM-iE(k)t/h\; (5) 

V Z7T J -co 



* out (x,t) = -= [f tr (k,t) + f ref (k,t)]e* kx dk, fout(k,t) = fUk,t) + f r e J(k,t); (6) 

V J-oo 



fUk,t) = ^T(k)A m (k)exp[z(J(k) - kd - E(k)t/h)}; (7) 



fZl(k,t) = ^R{k)A in {-k) exp[-i(J(fc) - F(k) - | + 2ka + E{k)t/K)] (8) 

where E(k) = h 2 k 2 /2m. 

For a completed scattering we have 

tyfuii(x,t) « Vi n (x,t) when t = 0, 

tyfuii(x,t) = ty out (x,t) when t -> oo. 

It is obvious that the larger is the distance a, the more correct is the approximation for 
f fu ii(x,t) at t = 0. 

For particles starting (on the average) from the origin, we have 

fok 

< x > in = 1 9 

m 

(hereinafter, for any Hermitian operator Q 

A _ < fin\Q\fin > 
< V >in~ , | , , 

^ Jin\Jin ^ 
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similar notations are used below for the transmitted and reflected wave packets). The 
averaging separately over the transmitted and reflected wave packets yields 

fit 

< x >L= - < k >L - < J'(k) >L +d; (fO) 



< x > r :j t = - < k >: e j t + < j\k) - F\k) +2« (n) 

m 

(hereinafter the prime denotes the derivative with respect to k). Exps. (9) — (11) yield 
the basis for defining the asymptotic tunneling times in the SWPA. 

2.2 Problems of the standard wave-packet analysis 

To display explicitly some shortcomings of the SWPA, let us derive again the SWPA's 
tunneling times. Let Z\ be a point to lie at some distance Li (Li ^> l and a — Li ^> / ) 
from the left boundary of the barrier, and Z 2 be a point to lie at some distance L 2 
(L 2 3> Iq) from its right boundary. Following [4], let us define the difference between the 
times of arrival of the CMs of the incident and transmitted packets at the points Z 1 and 
Z 2 , respectively (this time will be called below as the "transmission time"). Analogously, 
let the "reflection time" be the difference between the times of arrival of the CMs of the 
incident and reflected packets at the same point Z\. 
Thus, let t\ and t 2 be such instants of time that 

< x > in (h) = a - Li; < x > u out (t 2 ) = b + L 2 . (12) 

Then, considering (9) and (10), one can write the "transmission time" At tr (At tr — t 2 —t 1 ) 
for the given interval in the form 



m 

At tr = - 
n 



< k > tr 



out 



2 Li 

+ ha 

ko 



< k > r ut 



ko. 



(13) 



Similarly, for the reflected packet, let t[ and t' 2 be such instants of time that 

< x > in (t[) =< x > r ou { (t' 2 ) = a-L 1 . (14) 

From equations (9), (11) and (14) it follows that the "reflection time" At re f (At re / = 
t' 2 — t[) can be written as 



At 



ref 



m 

J 



<J'-F> > r ou {+U 



< -k >: ef 



out 



Li 

+ ^ + a 
ko 



. ^ ^out 



ko. 



(15) 
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Note that the expectation values of k for all three wave packets coincide only in the 
limit Iq — > oo (i.e., for particles with a well-defined momentum). In the general case these 
quantities are distinguished. For example, for a particle whose initial state is described 
by the Gaussian wave packet, we have 

A in (k) = Aexp(-l 2 (k- k ) 2 ), A = 
On can show that in this case 

< V > 

< k > tr = k + " ; (16) 

Let 

< k > tr = k + (Ak) tr , < ~k >ref= h + (Ak) re f, 

then relations (16) and (17) can be written in the form 

< T > m -(Ak) tr = -<R> m -(Ak) ref = <T .^ m . (18) 
Note that B! = -V . 

As is seen, quantities (13) and (15) cannot serve as characteristic times for a particle. 
Due to the last terms in these expressions the above times depend essentially on the 
initial distance between the wave packet and barrier, with L 1 being fixed. These terms 
are dominant for the sufficiently large distance a. Moreover, one of them must be negative. 
For example, for the transmitted wave packet it takes place in the case of the under-barrier 
tunneling through an opaque rectangular barrier. The numerical modelling of tunneling 
[1, 4, 5, 15] shows in this case a premature appearance of the CM of the transmitted packet 
behind the barrier, what points to the lack of a causal link between the transmitted and 
incident wave packets (see [8]). 

As was shown in [1, 4], this effect disappears in the limiting case l — > oo. For example, 
in the case of Gaussian wave packets the fact that the last terms in (13) and (15) tend 
to zero when l — > oo, with the ratio lo/a being fixed, can be proved with help of Exps. 
(16) and (17) (note that the limit l — > oo with a fixed value of a is unacceptable in this 
analysis, because it contradicts the initial condition a ^> l for a completed scattering). 
Thus, at first glance, in the limit l — > oo the SWPA seems to provide correct characteristic 
times for a particle. However, as will be seen from our formalism, even in this case, the 
above times are poorly defined. The point is that the above definitions of tunneling times 
for transmission and reflection are based on the implicit assumption that particles of both 
the subensembles start, on the average, from the origin as those of a whole quantum 
ensemble. As will be seen from the following, this is not the case. 
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Note, the fact that Exps. (13) and (15) cannot be applied to particles does not at all 
mean that they are erroneous. These expressions correctly describe the relative motion 
of the transmitted (or reflected) and incident wave packets. The principal shortcoming 
of the above approach is that it is meaningless to compare the motion of the transmitted 
(or reflected) wave packet with that of the incident one since they are related to different 
ensembles of particles and, as a consequence, there is no causal relationship between them. 
The right procedure of a separate timing of transmitted particles suggests the availability 
of such initial wave packet to evolve causally into the transmitted one. 

3 Separate description of transmission and reflection 
in the one-dimensional scattering problem 

For a long time the processes of transmission and reflection in a quantum scattering 
have been accepted to be inseparable, in principle. However, in this paper we show that, 
at least in the one-dimensional case, there are sound reasons to consider these processes 
separately. 

3.1 Wave function of a tunneling particle as a sum of wave func- 
tions to describe separately transmission and reflection 

According to quantum scattering theory, in one dimension, stationary-state wave functions 
for a particle impinging the barrier from the left (or from the right) possess one incoming 
and two outgoing waves. That is, in this case we deal, in fact, with one-source two-sinks 
scattering problems. 

Let for the problem at hand the amplitude of incoming wave, a in , be equal to unit, 
then the amplitudes of all four waves read as 



din 1) bout i a out j bi n (19) 

q q 

(note, B out = b out A in , A out = a out A in , B in = b in A in ). Let us also consider two auxiliary 
(two-sources one-sink) scattering tasks in which the amplitudes of incoming and outgoing 
waves are 



, re f 



iref 
u out 



q 



n re f 

a out 



o, 



uref 



P 



(20) 



and 



tr 



b tr 

u out 



o, 



L out 



1 
q 



b 



tr 



p* 



(21) 
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(the transfer matrix (3) is common for all three tasks). 

Note that in the first auxiliary task the only outgoing wave coincides with the reflected 
wave arising in the problem at hand (see (19)). And, in the second task, the only outgoing 
wave coincides with the transmitted wave in (19). It is evident that the sum of these two 
functions results just in that to describe the state of a particle in the initial tunneling 
problem. 

As is seen, the main peculiarity of the superposition of these two states is that due 
to interference the incoming waves, in the region x > b, disappear entirely (note that 
in the corresponding reverse motion they are outgoing waves). Figuratively speaking, 
interference reorients these waves into the region x < a. That is, in this superposition the 
probability fields of both sinks are radically reconstructed due to interference. Namely, 
they transform into fields with one outgoing and one incoming waves. 

Hereinafter, the wave function in which an incoming wave is associated with the re- 
flected wave of solution (19) will be refereed to as the reflection wave function (RWF). 
Similarly, the wave function in which an incoming wave is related to the transmitted wave 
of (19) will be refereed to as the transmission wave function (TWF). We postulate that, 
in the considered scattering problem, namely the nonstationary-state TWF and RWF 
describe, respectively, the transmission and reflection processes. As will be shown below, 
both the functions evolve in time with constant norms; at late times the TWF (RWF) 
coincides with the transmitted (reflected) wave packet. 

Thus, we see that the sum of wave functions (20) and (21) can be presented as that of 
the stationary-state RWF and TWF. Under the reconstruction of the probability fields, the 
squared amplitude of the incoming wave (in the region x < a) associated with reflection 
increases due to interference from the initial value \aZ f \ 2 (= R 2 ) (see (20)) to \a r ^\ 2 + 
Win\ 2 {= R 2 + TR = R) (in the RWF). In the case of transmission the corresponding 
quantity increases from the initial value \af n \ 2 (= T 2 ) (see (21)) to \af n \ 2 + \bf n \ 2 (= 
T 2 + TR = T) (in the TWF). 

Of course, the above postulate suggests the availability of a proper pair of solutions 
to the Schrodinger equation. The main thing which should be taken into account in 
finding these solutions is that the RWF describes the states of reflected particles only, 
and the TWF relates only to transmitted particles. As was said above, in both the cases, 
stationary solutions should contain one incoming and one outgoing wave. In this paper 
we show that such solutions do exist. 

3.2 Wave functions for one-dimensional transmission and reflec- 
tion 

So, let ty tr and ^ re f be the searched-for TWF and RWF, respectively. In line with 
subsection 3.1, their sum represents the wave function to describe, in the problem at 
hand, the sate of the whole ensemble of particles. Hence, from the mathematical point of 
view our task now is to find such solutions \l/ tr and \l/ re / to the Schrodinger equation that 
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for any t, 



*fuii(x, t) = V tr (x, t) + V ref (x, t) (22) 

where ^f u u(x,t) is the full wave function to describe all particles (see section 2). In the 
limit t — > oo 

*tr(a;,*) = < ut (x,0; *re/0M) = *) (23) 

where ^^(x, t) and ^^{(x, i) are the transmitted and reflected wave packets whose 
Fourier-transforms presented in (7) and (8). 

As is known, searching for the wave functions in the case of the time-independent 
potential V(x) is reduced to the solution of the corresponding stationary Schrodinger 
equation. For a given k, let us find firstly the functions *& re f(x; k) and ^tr( x ] k) for the 
spatial region x < a. In this region let 

* ref (x; k) = A m (k) (AZ f (k)e lkx + B r e J t (k)e- lkx ) (24) 

<M*; k) = A in (k)(AZ(k)e ikx + Bt t {k)e~^) (25) 

where Af n + A™ f = 1, % + B r J = b out . 

Since the RWF describes the state of reflected particles only the probability flux for 
ty re f(x] k) should be equal to zero, i.e., 

\AZ f \ 2 -\B r o e J\ 2 = 0. (26) 
In its turn, for ^tr( x ] k) we have 

\AZ\ 2 - \BHf = ^T(k) (27) 

lib 

(the probability flux for the full wave function tyf u u(x; k) and for ^(x; k) should be the 
same) . 

Taking into account that \l/ ir = ^ f u ii — ^ref let us now exclude \l/ tr from Eq. (27). As 
a result, we obtain for ty re f the equation 



Re (A^a\ n - BZ{b* out ) = 0. (28) 

The physical meaning of Eq. (28) is that the function ty re f(x), with zero probability flux, 
is such that the sum of the stationary-state RWF and any other stationary-state wave 
function with a nonzero probability flux does not change the value of the latter. 
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From condition (23) for ^ ref (x; k) it follows that B r Q e J t (k) = b out (k) = p*/q (see (19)). 
Then Eq. (28) yields Re(A r ^) = R, and Eq. (26) leads to \A™ f \ 2 = \B r e u {\ 2 = \p*/q\ 2 = R- 
Thus, AZ f = VR{VR±i\ / T) = v^Rexp(iA); A = ± arctan( ^T/i?). 

So, there are two solutions to satisfy the above requirements for ^/ re f(x; k), in the 
region x < a. Considering Exps. (4) for the elements q and p, we have 



^ref{x; k) = -2\fRA in sin (k(x - a) + -(\ - J + F - )e' 



<£(+) 



(29) 



where 



(±) 



\±[J-F-- + 2ka) 



Now we have to show that only one of these solutions describes the state of the 
subensemble of reflected particles. To select it, we have to study both the solutions in the 
region x > b where they can be written in the form 



* ref (x; k) = A m {k) (Kj t {ky kx + BZ f (k)e-^) 



(30) 



where 



Considering Exps. (4) as well as the equality exp(iA) = \f~R ± iVT, one can show that 

G = =p exp i(kb - ^ J + F + ^ - A^ ; 
here the signs (=f) correspond to those in the expression for A. Then, for x > b, we have 



^ re/ (x; k) = Tl\fRA in sin \k(x -b) + ]-(j + F+^--\ 



(31) 



For the following it is convenient to go over to the variable x'\ x = x mi d + x' where 
Xmid = (a + b)/2. Then we have, for x' < —d/2, 



m ref (x') = -2\fRA in sin [hkd + A - J - ^) + ^ + kx' 



for x' > d/2 



^ref(x') = ±2\fRA in sin [^(kd + A - J - |) - ^ - kx' 



i( t>(+) 



i<P(+) 



From these expressions it follows that for any point x' = x (x < —d/2) we have 



^ref(^o) = -2\[RA in sin -{kd + A - J h F) + kx ( 



^(+) 



(32) 
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^ref(-xo) = ±2v r RA in sm \-(kd + \- J- ^- + F) + kx - F]e^(+). (33) 

Let us consider the case of symmetric potential barriers: V(x') = V(-x'). For such 
barriers the phase F is equal to either or ir. Then, as is seen from Exps. (32) and (33), 
one of the above two stationary solutions \l/ re /(x'; k) is odd in the out-of-barrier region, 
but another function is even. Namely, when F = the upper sign in (33) corresponds to 
the odd function, the lower gives the even solution. On the contrary, when F = 7r the 
second root A leads to the odd function ^ re f(x'; k). 

It is evident that in the case of symmetric barriers both the functions keep their "out- 
of-barrier" symmetry in the barrier region as well. Thus, the odd solution \l/ re j(x'; k) is 
equal to zero at the point x' = 0. Of importance is the fact that this property takes 
place for all values of k. In this case the probability flux, for any nonstationary-state 
wave function formed only from the odd (or even) stationary solutions ^/ re f(x'; k), should 
be equal to zero at the barrier's midpoint. This means that for particles impinging a 
symmetric barrier from the left they are reflected by the barrier without penetration into 
the region x' > 0. In its turn, this means that the searched-for stationary-state RWF 
should be zero in the region x' > 0, but in the region x' < it must be equal to the 
odd function ^/ re f(x'; k). In this case the corresponding probability density is everywhere 
continuous, including the point x' = 0, and the probability flux is everywhere equal to 
zero. 

Of importance is the fact that the above property of reflection admits, in principle, 
experimental checking. Indeed, since reflected particles does not penetrate into the region 
x > x rni( i of the symmetric barrier, the switching on an infinitesimal magnetic field in this 
region must not influence the spin of these particles. For checking this property, one can 
use the experimental scheme presented in [19]. 

As regards the searched-for TWF, \l/* r (x; k), it can be found now from the expression 
^tr{%', k) = tyf u u(x;k) — ^/ re f(x;k). This function is everywhere continuous, and the 
corresponding probability flux is everywhere constant (we have to stress once more that 
this quantity has no discontinuity at the point x = x mic i, though the first derivative of 
^tr(x; k) is discontinuous at this point). Thus, as in the case of the RWF, wave packets 
formed from the stationary-state TWF should evolve in time with a constant norm. 

As is seen from Exps. (32) and (33), for asymmetric potential barriers, both the 
solutions ^ re f(x'; k) are neither even nor odd functions. Nevertheless, it is evident that 
for any given value of k one of these solutions has opposite signs at the barrier's boundaries. 
This means that, for any k, there is at least one point in the barrier region, at which this 
function is equal to zero. However, unlike the case of symmetric barriers, the location of 
such a point depends on k. Therefore the behavior of the nonstationary-state RWF in the 
barrier region is more complicated for asymmetric barriers. Now the most right turning 
point for reflected particles lies, as in the case of symmetric barriers, in the barrier region, 
but this point does not coincide in the general case with the midpoint of this region. 



13 



To illustrate the temporal behavior of all the three wave functions, i.e., ^ full-, ^tr and 
ty re fi we have considered the case of rectangular barriers. In this case, the stationary-state 
wave function ty re f(x] k), for a < x < x mi d, reads as 

^ ref = 2\ r RA in e i ^+'>[cos(ka + sinh(/trf/2) 
k 

sin(A;a + 0r_)) cosh(/tc?/2)] sinh(/t(x — x mic i)) (34) 

K 



where k = w2m(t> — E)/h (the below-barrier case); and 



^ ref = -2\fRA in e i4) w[cos{ka + ( _ } ) sin(/td/2) 
A; 

+ — sin(A;a + </>(-)) cos(/tc?/2)] sin(/t(a; — x m idj) (35) 



where n = ^2m{E — v )/h (the above-barrier case). In both cases ^/ re f(x; k) = for 

We have calculated the spatial dependence of the probability densities \^f u u(x, t)\ 2 
(dashed line), \^f t r(x, t) | 2 (open circles) and \^ re f(x, t) | 2 (solid line) for the rectangular 
barrier (Vo = 0.3eK, a = 500nm, 6 = 505nm) and well (T^o = — 0.3eV, a = 500nm, 
6 = 505nm). Figures 1 (t = 0), 2 (t = OAps) and 3 (t = 0.42ps) display results for the 
barrier, and figures 4 (t — 0), 5 (t — OAps) and 6 (t — 0A3ps) display results for the 
well. In both the cases, the function tyf u u(x, 0) represents the Gaussian wave packet with 
l = 7.5nm; the average kinetic energy is equal to 0.25eV, both for the barrier and well. 
Besides, in both cases, the particle's mass is 0.067m e where m e is the mass of an electron. 

As is seen from figures 1 and 4, the average starting points for the RWF and TWF 
differ from that for ^ j u ii- The main peculiarity of the transmitting wave packet is that 
it is slightly compressed in the region of the barrier, and stretched in the region of the 
well. Figure 7 shows that, at the stage of the scattering event it = OAps; see also figure 
2), the probability to find a transmitting particle in the barrier region is larger than in 
the neighborhood of the barrier. This means that in the momentum space this packet 
becomes wider when the ensemble of particles enters the barrier region. For the well (see 
figure 8) there is an opposite tendency. Note that for the barrier < T >j„~ 0.149. For 
the well < T > in ^ 0.863. 



3.3 Connection of the wave functions for reflection and trans- 
mission with the eigenvectors of the scattering matrix 

Of importance is the fact that there are other two settings of the tunneling problem for 
the given potential V(x) when the subensemble's states described by the RWF and TWF 
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arise explicitly. Indeed, let us find such solutions to the Schrodinger equation, for a given 
potential V(x), for which 

(t:H(t) <*» 

where S is a constant. This means that the amplitudes of incoming waves should obey 
the characteristic equation 





s = (fk (37) 

where S is the scattering matrix. 

It is easy to show that the solutions of this equation can be written in the form 



S = 



l + ifJ.\p\ f a in \ _ ( i/ip/\p\ 

\ 1 




C (M) 



where fj, — ±1; C(+) and C(_) are arbitrary constants. 

Now let us find such values of C(+) and C(_) at which 6 out = p*/q. It easy to show that 
all four amplitudes read, in this 

l+tfi\p\ q 

yM s m.l. t bln = ^?=jLvR(VR + inVT). (38) 
q p q 1 + iH\p\ i>lAp\ 

One of two solutions with these amplitudes is evident to coincide, for x < a, with the 
RWF found in subsection 3.2. This means that in the case of symmetric potential barriers 
this function, like the RWF, is equal to zero at the midpoint of the barrier region, for 
any value of k. In this two-sources scattering problem, both the incident wave packets 
does not cross the above point. In fact, we deal here with the ideal bilateral reflection of 
particles from the midpoint of the barrier region, which is described by the sum of two 
the RWFs. 

In a similar way, for the same eigenvalue of the scattering matrix, one can find such 
values of C( + ) and C(_) at which a out = 1/q: 



Vr(Vr - i/iVR); b out = • -; 



1 + ifi\p\ p q 

a ut = -; b m = - 1 , , = • VT(Vt - iptVR). (39) 

q p 1 + tfi\p\ p 
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As is seen, the stationary-state TWF appears explicitly in the solution with amplitudes 
(39). In the case of symmetric potential barriers this solution is evident to represent 
a sum of two continuous wave functions whose probability fluxes are continuous too. 
For one of them the amplitudes of incoming and outgoing waves are, respectively, a,i n 
(= \/T(\fT — i/iy/R)) and a out (= l/q). For another function these amplitudes are, 
respectively, b in (= -(ifJ,\p\/p) ■ VT(VT - ifxy/R)) and b out (= -(ifJ,\p\/p) ■ (l/q)). The 
first (second) function is just the TWF to describe the ideal transmission of particles 
impinging the barrier from the left (right). The corresponding nonstationary-state wave 
functions are evident to evolve in time with a constant norm. 

So, each of the above "two-sources" wave functions generated by eigenvectors of the 
scattering matrix represent a sum of two causally evolved "one-source" wave functions. 
One of them describes the state of a particle impinging the barrier from the left. Another 
function relates to particles moving to the right of the barrier. In the case of (38) both 
the one-source wave packets are ideally reflected by the barrier. And, in the case of (39) 
both one-source wave packets are ideally transmitted by it. These two auxiliary tunneling 
problems give us the basis to verify the formalism presented in this paper. 

Note also that the stationary-state RWF and TWF, for the problem at hand, should 
correspond to the same value of /i, i.e., to the same eigenvalue of the scattering matrix. 
As regards another eigenvalue, in the case of reflection it generates the even function 
which does not fit as a RWF (see subsection 3.2). That is, only one of the eigenvalues 
of the scattering matrix is associated with the RWF and TWF of the scattering problem 
considered. 



4 Exact and asymptotic tunneling times for trans- 
mission and reflection 

4.1 Exact tunneling times 

So, we have found two causally evolved wave packets to describe the subensembles 
of transmitted and reflected particles in the considered tunneling problem, at all stages 
of the scattering process. As is shown, the motion of these packets can be, in principle, 
observed experimentally. It is evident that the given formalism may serve as the basis 
to solve the tunneling time problem, since now one can follow the CMs of wave packets, 
which describe separately reflection and transmission, at all instants of time. 
Let tf and t\ be such instants of time that 



< ^ r (x,tf)|^ r (Mf) > 
< m tr (x^)\x\m tr (x,tf) > 
<V tr (x,t t {)\y tr (x,t t {)> 



a-L x - (40) 
b + L 2 , (41) 
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where ^/ tr (x,t) is the subensemble's wave function found above for transmission. Then, 
one can define the transmission time At tr (L 1: L 2 ) as the difference t 2 (L 2 ) — ^(-^i) wnere 
t l {(Li) is the smallest root of Eq. (40), and t 2 r (L 2 ) is the largest root of Eq. (41). 
Similarly, for reflection, let ti(Li) and i£ (£1) be such instants of time t that 

< V ref (x,t)\x\V ref (x,t) > 

Z ; TTT ; s — Oi — Li\. l^zi 

<V ref (x,t)\y ref (x,t)> 

Then the reflection time Ai re /(L 1 ) can be defined as At re f(Li) = f 2 f — t[ ef where t™* {Li) 
is the smallest root, and f 2 (L 2 ) is the largest root of Eq. (42) (of course, if they exist). 

It is important to emphasize that, due to conserving the number of particles in both 
the subensembles, both these quantities are non-negative for any distances Li and L 2 . 
Both the definitions are valid, in particular, when L\ = and L 2 = 0. In this case the 
quantities At ir (0,0) and At re /(0) yield, respectively, exact transmission and reflection 
times for the barrier region. Of course, one has to take into account that in the case of 
reflection the CM of the wave packet may turn back without entering the barrier region. 



4.2 Asymptotic tunneling times 

It is evident that in the general case the above average quantities can be calculated only 
numerically. At the sane time, for sufficiently large values of L\ and L 2 , one can obtain 
the tunneling times At tr (Li, L 2 ) and At tr (Li, L 2 ) in more explicit form. Indeed, in this 
case, instead of the exact subensemble's wave functions, we can use the corresponding 
in asymptotes derived in A;- representation. Indeed, now the "full" in asymptote, like the 
corresponding out asymptote, represents the sum of two wave packets: 

/£(M) = ^r{k)A m (k)exp[t(A(k)-a^-E(k)t/h)}; (43) 

f[: f (k,t) = jR(k)A in {k) exp[*(A(A;) - E(k)t/h)]; (44) 

a = 1 if A > 0; otherwise a = —1. Here the function A(k) coincides, for a given k, with 
one of the functions, X(k) or — X(k), for which ty re f(x] k) is an odd function (see above). 
One can easily show that for both the roots 

IT'I 

|A'(fc)| = -^=L. 

A simple analysis in the /c-representation shows that well before the scattering event 
the average kinetic energy of particles in both subensembles (with the average wave num- 
bers < k >f n and < k > r in ^) is equal to that for large times: 

< k -> tr —< k -> tr < k -> re f— — <r k > re/ 
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Besides, at early times 



fit 

< * >Z= ~<k>Z-< A'(k) >Z; (45) 



< £ >::/= -<k>z f -< M(k) >z f (46) 

As it follows from Exps. (45) and (46), the average starting points x% arl and x r s tLt, 
for the subensembles of transmitted and reflected particles, respectively, differ from that 
for all particles: 



x tr 

start 



< A'(fc) >t, x r s tt = -< A'(fc) >Z f • (47) 



The implicit assumption made in the SWPA that incident, as well as transmitted and 
reflected particles start, on the average, from the same point does not agree with this 
result. By our approach, this is the main reason why the asymptotic transmission and 
reflection times obtained in the SWPA should be considered as ill-defined quantities, for 
any wave packets. 

Let us take into account Exps. (45), (46) and again analyze the motion of a particle 
in the above spatial interval covering the barrier region. In particular, let us calculate the 
transmission time, r tr , spent (on the average) by a particle in the interval [Zi,Z 2 ]. It is 
evident that the above equations for the arrival times t\ and t 2 r , which correspond the 
extreme points Z\ and Z 2 , respectively, read now as 

< x >t n (if) = o - Li; < x >* t (#) =b + L 2 . 
Considering (45) and (10), we obtain from here that now the transmission time is 

T tr(Li, L 2 ) EE 4 - t\ r = h< m k> tr n (< ^ >out " < A' >f n +Z* + L 2 ) . (48) 

Similarly, for the reflection time, r re f(Li) (r re f = t 2 e * — f^), we have 

< x >Z f (t? f ) = a-L 1 , <x > r e u { (f 2 s ) =a-L 1 . 
Considering (46) and (11), one can easily show that 

r re/ (Li) = f 2 ef - t[ ef = ^ (< J' - F> > r ; u { - < A' >Z f +2Lx) . (49) 

The inputs (t™ = r tr (0, 0)) and r™f {t 1 ^ = r tr (0, 0)) will be named below as the 
asymptotic transmission and reflection times for the barrier region, respectively: 
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< s = n ™k>& ( K J ' >L ~ < A ' >Tn ) ' (50) 

C/ = r ™ re/ (< ^ - ^ > r :( - < a' >i: / ) (si) 

« < « >m 

Here the word "asymptotic" points to the fact that these quantities were obtained with 
making use of the in and out asymptotes for the subensembles investigated. Unlike the 
exact tunneling times the asymptotic times may be negative by value. 
The corresponding lengths df^ and cTjff, 

€ff =< J' >out - < A' >Tn, €/f =< J' - F' >Z f t - < A' >Z f , (52) 

can be treated as the effective widths of the barrier for transmission and reflection, re- 
spectively. 

4.3 Average starting points and asymptotic tunneling times for 
rectangular potential barriers and ^-potentials 

Let us consider the case of a rectangular barrier (or well) of height Vq and obtain explicit 
expressions for d e ff(k) (now, both for transmission and reflection, d e ff(k) = J'(k) —A'(k) 
since F'(k) = 0) which can be treated as the effective width of the barrier for a particle 
with a given k. Besides, we will obtain the corresponding expressions for the coordinate, 
Xstart(k), of the average staring point for this particle: x sta rt(k) = —A'(k). It is evident 
that in terms of d e ff the above asymptotic times for a particle with the well-defined 
average momentum fc read as 

T as = as = md eff (ko) 

tr ref hk ' 

Using the expressions for the real tunneling parameters J and T (see [36, 39]), one 
can show that, for the below-barrier case (E < V ), 
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k 2 + nl sinh 2 (nd/2) [k,q smh(nd) - k 2 Kd] 



k Ak 2 K 2 + Kq sinh (red) 

, . (k 2 — k 2 ) smh(Kd) + k 2 nd cosh(/trf) 

XstartW = -2- 

where k = y / 2m(Vb — E)/fi 2 ; for the above-barrier case (E > Vq) - 



4 

d eff(k) = - 



k 2 - (3k 2 sin 2 («d/2)] \k 2 Kd - (3k 2 , sin(/cd)] 



k 4k 2 k 2 + Kq sin (red) 
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(k 2 + k 2 ) sm(K,d) — k 2 nd cos(nd) 

X sta rt{k) = -2(3- 

where k = pm(E - V )/h 2 ; (5 = 1 if V > 0, otherwise, (3 = — 1. In both the cases 
k = ^2m\V \/h 2 . 

It is important to stress that, in the limit k — > oo, d e ff d and x s t ar t{k) — > 0. This 
property guarantees that for infinitely narrow in x-space wave packets the average starting 
points for both subensembles will coincide with that for all particles. It is important also 
that for wells the values of d e ff and, as a consequence, the corresponding tunneling times 
are negative, in the limit k — > 0, when sm(n d) < 0. 

Note that for sufficiently narrow barriers and wells, namely when /td <C 1, we have 
d e ff ~ d. That is, particles tunnel, on the average, classically through such barriers. For 
the starting point we have 

K 2 K 2 
Xstart(k) W -^d, Xstart(k) « ~l 3 ^. d 

for E < V and E > V , respectively. 

For wide barriers and wells, when nd ^> 1, we have d e ff ~ 2/k and x start (k) ~ 0, for 
E < V ; and 

,,9, k 2 - 8 sin 2 (K,d/ 2) .,. 28K 2 l k 2 dcos(Kd) 

d ef f W 4Frf • ; " , V 7 ( X star t{k) W ■ , 2 ", 4 ■ 2 , ' , 

for ,B > Vb. 

It is interesting to note that for the 5-potential V(x) = W5(x— a) d e ff = 0. This means 
that, contrary to the phase tunneling time, the tunneling times defined here equal to zero 
for this potential. As regards the starting point x star t(k) in the case of the (^-potential, 
we have 

2mh 2 W 

Xstart{k) - - h i k2 + m2w2 - 

Thus, we see that, for example, in the case of 5-wells (W < 0) particles in each subensem- 
ble start, on the average, with an advance in comparison with those of the whole quantum 
ensemble. 



5 Conclusion 

A separate description of transmission and reflection is commonly accepted to contradict 
the principles of quantum mechanics. However, in this paper we argue that this is not the 
case, at least in the one-dimensional one-particle scattering problem. We show that the 
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wave function to describe, in this problem, the state of the whole ensemble of identically 
prepared particles can be uniquely presented as the sum of two functions (named here 
as the TWF and RWF) to obey the Schrodinger equation. In the case of nonstationary 
case, these functions are mutually orthogonal. At late times the TWF coincides with the 
transmitted wave packet, and the RWF approaches the reflected one. We postulate that 
namely the TWF and RWF are the wave functions to describe, respectively, transmission 
and reflection in the considered scattering process. 

There is also a widely accepted viewpoint (see, for example, page 106 in [40] and page 
17 in [41]) that all solutions to the stationary one-dimensional Schrodinger equation, for 
a finite potential, must be everywhere continuous together with their spatial derivatives; 
otherwise, the points where this requirement is violated contain allegedly sinks or sources 
of particles. However, in this paper we show that the above requirement for "physical" 
solutions to the Schrodinger equation is, in reality, excessive. The main peculiarity of the 
presented stationary-state wave functions for transmission and reflection is that there is a 
point in the barrier region where these everywhere continuous functions have discontinuous 
first spatial derivatives (in the case of symmetric potential barriers this takes place at the 
midpoint of the barriers). Nevertheless, for each subensemble, this point contains neither 
sink nor source of particles: both for transmission and reflection, the probability current 
density for each stationary-state wave function is constant on the spatial axis, and the 
norm of wave packets formed from these functions is constant in time. 

We show that, in the case of a symmetric potential barrier, reflected particles impinging 
the barrier from the left do not penetrate into the spatial domain lying to the right of 
the midpoint of the barrier region. This means, in particular, that the switching on an 
infinitesimal magnetic field in this domain must not influence the spin of these particles. 

Besides, for the given potential we formulate two scattering problem in which the 
RWF and TWF arise separately and, as a consequence, there is another possibility to 
check experimentally our approach. In both the scattering problems the amplitudes of 
incoming waves form the eigenvectors of the scattering matrix for the given potential. 

On the basis of the above formalism we define average (exact and asymptotic) trans- 
mission and reflection times. The exact tunneling times are always non-negative. In the 
case of rectangular barriers and 5-potentials, for both the subensembles, we derive explicit 
expressions for the asymptotic tunneling times and for the average starting points. These 
times differ essentially from those arising in the SWPA. 
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Figure captions 

Fig. 1 The x-dependence of \^f u a{x, t)\ 2 (dashed line) which represents the Gaussian 
wave packet with l pa 7.5nm and the average kinetic particle's energy 0.25eV, as well 
as \^ tr (x, t)\ 2 (open circles) and \^ re f(x, t)\ 2 (solid line) for the rectangular barrier (V = 
0.3eV, a = 500nm, b = 505nm); t = 0. 

Fig. 2 The same as in Fig. 1, but t = OAps. 

Fig. 3 The same as in Fig. 1, but t = 0A2ps. 

Fig. 4 The x-dependence of \^f u u(x, t)\ 2 (dashed line) which represents the Gaussian 
wave packet with l = 7.5nm and the average kinetic particle's energy 0.25eV, as well 
as \^tr(x, t)\ 2 (open circles) and \^ re f(x, t)\ 2 (solid line) for the rectangular well (V = 
— 0.3eV, a = 500nm, b = 505nm); t = 0. 

Fig. 5 The same as in Fig. 4, but t = OAps. 

Fig. 6 The same as in Fig. 4, but t = 0A3ps. 

Fig. 7 The same functions for the barrier region; parameters are the same as for Fig. 2. 
Fig. 8 The same functions for the barrier region; parameters are the same as for Fig. 5. 
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